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Abstract 


Scintillations are observed during occultations of both stars and spacecraft 
by planetary atmospheres. Existing treatments of these scintillations ignore a 
major effect unique to occultations: the severe flattening of the Fresnel zone 

or source image by defocusing. Other large effects, due to ’’saturation" of the 
scintillation, have also been ignored. The deeper portions of atmospheric 
temperature and density profiles inferred from occultation data are seriously in 
error if other planets' atmospheres are as turbulent as our own. Thus, profiles 
obtained from entry probes (e.g,, the Soviet Venera series) are probably more 
accurate than those from radio occultation (Mariner 5 and 10) data. Scintillation 
greatly reduces the information obtainable from occultation observations; much of 
the detail in published profiles is probably spurious. This paper gives an 
approximately - correct theoretical treatment that is a substantial improvement 
over published theories, and shows how a more accurate theory could be constructed. 
Some methods for a more accurate determination of atmospheric structure are proposed. 
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I. Introduction 


In recent years, a large literature has grown up around the inference of 
atmospheric structure (e.g., temperature and pressure profiles) from observations 
of the intensity and/or phase of electromagnetic waves received from a point 
source, such as a star or a spacecraft radio transmitter, as it is occulted by a 
planetary atmosphere. The occultations of Regulus by Venus in 1959, of BD - 
1704388 by Neptune in 1968, of 8 Scorpii by Jupiter in 1971, and of Mariners 5 
and 10 by Venus in 1967 and 1974 respectively, as well as numerous occultations 
of other Mariner spacecraft by Mars*, have been used to infer planetary atmospheric 
profiles. Many events were observed at several stations, and some observations 
have been re-discussed several times, so the total number of publications is very 
large; I shall not attempt to review this extensive literature. 

The theoretical interpretations of such data go back at least to Pannekoek 
(1903), who wished to explain the phenomena of the occultation of a star by Jupiter 
on Sept. 19 of that year. He commented that the light curves of future occultations 
of stars by planets should be measured accurately: "This is of importance, because 
It can give a determination of the temperature distribution in the outermost 
layers of the planetary atmosphere, or more precisely, the way in which the hori- 
zontal refraction varies with height." The problem was re-analysed (in French 
instead of in German) by Fabry (1929), who also considered the effects of both 
diffraction and an extended source, and vho drew the apparent shape of the occulted 
object at various stages. A third derivation of the laws of occultations was 
published in English, with a third independent scheme of notation, by Baum and 
Code (1953) , who had observed the occultation of a Arietis by Jupiter in the pre- 
vious year. Finally, the subject was reviewed by Link (1969), using a fourth set 
of symbols. Since then still other notations have been used by various authors. 
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This lack of uniformity is a problem, complicated here by the need tc discuss 
scintillation as well (see Table I). 

All the above authors assume a simple exponential atmosphere with a fixed 
scale height. More recently, it has become fashionable to reduce occultatlon 
curves by a numerical inversion process, which assumes a continuous, one-to-one 
mapping of rays between the successive layers of the planet's atmosphere and the 
successive moments of observation (or points along the plane of the telescope 
objective.) As will be shown below, this assumption is probably not justified, 
at least for the deeper stages of the occultation, and therefore the inversion 
results contain systematic errors. 

The first hint that turbulence might play a part in the observed phenomena 
was a comment by Baum and Code (1953) : "A sequence of photographs taken by 

Pettit and Richardson with the 100-inch telescope indicates that the fading is 
far from smooth and in fact suggests that the principal humps in [the light curve] 
may partly represent real features in the occultation - curve. Random variations 
in 0 caused by turbulence in Jupiter’s stratosphere would tend to jumble neigh- 
boring rays in such a way that the observed flux . . . would deviate randomly from 
an ideal occultation-curve. Percentage variations in 0 of the same order as 
those suffered by objects near the earth’s horizon v.’ould become significant for 
0 greater than the angular subtense of Jupiter’s scale height ... which occurs 
when <J> < 0.5 <£ 0 . Thus, large fluctuations in the tail of the observed light- 
curve are rather to be expected.” 

Osawa et al. (1968) remarked that their light-curves "have large fluctua- 
tions prevailing the main ptuiOc of occultation. Small fluctuations lasting for 
about two minutes after immert- . o before emersion are also remarkable fea- 
tures of these light-curves. Random variations in refraction angle caused by 
turbulence in Keptune’s upper atmosphere must be the cause of these fluctuations.” 
They did not, however, take the important step of showing that realistic amounts 
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of turbulence cculd account for their observations, as Baum and Code (1953) had 
done. 

Freeman and Lynga (1969, 1970) discussed their observations of the same 
event. They seem to have been the first workers to call the observed asymmet- 
rical light fluctuations ’'spikes*' , and also to have attributed them to layers of 
lower-than-average refractivity gradient in the planetary atmosphere. Since 
then, this explanation has become widely accepted. According to this model, 
"spike-"'" in the light curve correspond to layers with a more nearly adiabatic 
lapse rate, and the minima between them to more isothermal or stably - stratified 
regions. 

This view was endorsed, for example, by Kovalevsky and Link (1969), the 
first workers to use the numerical-inversion technique instead of curve-fitting. 
They claimed, without any quantitative demonstration, that the fluctuations 
observed at two stations several hundred kilometers apart showed "a good 
correlation." 

Unfortunately, the human eye tends to "see" patterns in random distributions; 
the subjective impression of "a good correlation" may not be justified by the 
actual data. For example. Rages et al. (1974) concluded from the same data as 
Kovalevsky and Link (1969) that "the. correlation is not striking." Furthermore, 
in a quantitative statistical discussion. of spikes in the 6 Sco occultation, 
Veverka et al. (1972) concluded, "We find no evidence in our data that the pattern 
of fluctuations is planet-wide in the case of Jupiter, and it is therefore im- 
proper to speak of global layering." 

Finally, although Veverka et al. (1974a) assume that the "spikes" are due to 
"layering," and claim that "the agreement among the three channels in Fig. 10 
[their temperature-ref ractivity profiles] is excellent," closer inspection of 
this figure and their Fig. 13 (which shows refractivity scale height as a function 
of depth) shows that the fine structures of the atmospheric profiles deduced by 
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inversion from observations at one station in different wavelengths are not at 
all in good agreement. While the same general features appear in these profiles, 
th.?2 do not occur at the same depth in the atmosphere; in many cases a locrl 
maxitaim on one curve occurs at the same atmospheric level as a minimum on another 
curve. This discrepancy is explained by the scintillation model below. 

A similar problem is shown by the very disparate profiles deduced by Rages 
et al. (1974) from three different observations of the Neptune event. Their pro- 
files differ by several tens of degrees at point? a few hundred kilometers apart 
on Neptune. It is difficult to see how a planet so far from the Sun can develop 
such large horizontal temperature gradients, which (if real) would cause very 
strong thermal winds . 

Other problems with temperature pt -files deduced from occultations are 
reported by Veverka and Washerman (1974) and Kliorc et al. (1974). Both the 
Regulus and the Pioneer 10 oceultation data lead to temperature profiles that are 
systematically too high, by roughly 100°K. The discrepancies are so striking 
that Veverka and Vasserman (1974) soy, "There is no way in which the Regulus 
oceultation light curve can be reconciled with reality," and Kilore et al. say 
u lhis temperature distribution in the upper levels of the atmosphei'e of Jupiter 
is completely different from any published model of the atmosphere," 

These statements cast doubt on the validity of numerical-inversion results. 
The spurious results are, 1 believe, due to scintillation effects. 

II. Scintillations: Theory 

Although it has been common practice to interpret the fluctuations in signal 
strength at radio frequencies as scintillations, using the theory for isotropic 
Kolmogorov turbulence, this has not yet been applied to the optical observations. 
Furthermore* the theory' that has been used (cf. Woo and Ishimaru, 1974) is appro- 
priate to a medium without differential refraction, in which t-;-> Fresnel zone 


is circular. Finally, only weak-scintillation theory has been used, although 
many of the observations show very strong modulation. Thus the first need is to 
adopt scintillation theory to suit the special circumstances of planetary 
occultations. 

1. Effects of Differential Refraction. 

The theoretical treatments used up to now have implicitly assumed a circular 
Fresnel zone. But, as Fjeldbo and Eshleman (1969) point out, "Refraction has the 
effect of changing these zones from circular to elliptic. In the atmospheric and 
ionospheric regions that produce defocusing, the curvature of the wavefronts is 
increased in such a way as to squash the Fresnel zones into ellipses with minor 
axes in the vertical direction." This is strictly true only if the Fresnel zone 
is much smaller than a scale height; if it is comparable, the Fresnel zone (or 
the disk of an occulted star, whichever is larger) is more flattened on the 
lower side. The shapes that result are illustrated in Fabry (1929), Link (1969), 
or the form of the setting Sun - a most familiar example of an atmospheric 
occultatiun with differential refraction (Young, 1974). In our present examples, 
the free-space Fresnel zone sizes are roughly a kilometer across, and the atmos- 
pheric scale heights are typically ten times larger. Thus the elliptical approxi- 
mation is fairly good, but not exact. 

From the basic physical law that the power received at a point is just what 
flows through the first Fresnel zone, one sees that the reduction in area of this 
zone (or the axial ratio, in the elliptical approximation) is just the factor «*> 
by which the intensity has been reduced. For an extended source (stellar disk) 
the corresponding result follows from the fact that a lossless optical system 
(the occulting atmosphere) cannot change the source's apparent surface brightness; 
so the defocusing produces an image-area reduction equal to the reduction in 
received intensity. Thus, when the received signal is cut in half, the region 
of the atmosphere probed is half as large vertically as horizontally (see Fig. 1) . 


- 8 - 


This reduction in size has a profound effect on the scintillations, which are 
mainly due to the smallest elements of turbulence that arc not averaged out over 
either the source or the Fresnel-zonc area. 

To estimate this effect quantitatively, an heuristic approach used previ- 
ously (Young 1969, 1970) is useful. For detailed accounts of scintillation 
theory, see Tatarskii (1961, 1967), Lee (19t>9), or Lee and Harp (1969). We can 
regard the effect of either Fresnel-zone or source size as an averaging aperture 
in the turbulent medium, and then use the ray-optical model (Reiger, 1963; 

Young, 1967, 1969) of scintillation. Then the mean-square modulation (or "scin- 
tillation power," P) of the shadow pattern observed in the x-y plane Is propor- 
tional to 

p - f* dk x dk / dz • C>) k 1/3 z 2 F (k, z, X, ...), (1) 

2 2 U 

where k - (k + k ) is the spatial frequency of a Fourier component of the 
x y 

shadow pattern, z is the distance from a slab of turbulence to the observer, 

C (z) is a measure of the strength of the turbulence at distance z, X is the 
n 

wavelength of the light, and F (...) is a spatial frequency filter function 
corresponding to source or Fresnel zone size. The effect of this filter is to 

cut off the infinite integrations in k-space at finite values (Young, 1970). 

1/3 2 
The factor k ' is due to the Kolmogorov spectrum of turbulence, and z is due 

to a "lever-arm" effect (Young, 1969) that represents the gradual conversion of 
phase fluctuations in the wavefront to amplitude fluctuations. 

We can readily obtain the temporal spectrum of the scintillations (Reiger, 

1963; Young, 1967, 1969; Lee and Harp, 1969) by choosing the x-axis in the direc- 
tion of motion of the line of sight through the turbulent medium, omitting the 
integration over k , and transforming from spatial frequency k in the direction 
of motion to temporal frequency f: 

f e v l k x * 

2u 


( 2 ) 
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where is the speed of motion. The integrals over z and k then give P (k ) 

y 

or P (f), the spectral power density. For a point source at infinity, P is the 
"diffraction filter" function (Young, 1969, 1970) 


diff 


sin*' (Azk 2 /4x) 


( 3 ) 


(Azk 2 /4x) 2 


~k 

which cuts off at k : (Az) . The high-frequency behavior of the spectral power 

e ~8/3. 


density l’(f) is then readily found to be asymptotically proportional to f ; 
a simple expansion for calculating the spectrum from a single slab of turbulence 
is given by Young (1971) . For a uniform disk source of angular radius r, the 


averaging region at z is a circle of radius rz, and the filter is (Young, 1969) 

,2 


F di6k « [2J 1 (rzk)/(rzk)]‘ 


(4) 


-1 


where is the usual Bessel function; this cuts off at k u (rz) . Lee (1974) 


has discussed the effects of other interesting source shapes. 

For a region of arbitrary shape, the filter F has the analytical form of 
the intensity distribution in its Fraunhofer diffraction pattern. Thus, for an elliptic 
averaging region with semimajor axis a and semiminor axis $a, where £ 1, we have 

F cU = [2J l (ak*)/(al:*)] 2 , (5) 

where 


k* * t(k /*) 2 + k 2 ) ** , (6) 

a y 

assuming the major axis in the y-direction. That is, compressing the averaging 
region in the x-direction by a factor ^>1 has effectively extended the cut- 

off of the filter function by this same factor along the k^ axis. In other words, 
when we substatute Equations (5) and (6) into Equation (1), the k^ integral cuts 
off at ky = a \ and the k^ integral cuts off at k^ = (a<j) = i>/a. In the 
occult.ation light curve, $ is the ratio of the partially - occulted intensity to 
the unocculted intensity, and is typically much smaller than unity when scintilla- 
tions are observed. An atmospheric averaging region that is very different in the 
vertical and horizontal directions results from atmospheric dispersion on Earth 
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(Young, 1970), and the same treatment can be used here; the chief difference is 
that the sense of the asymmetr’ is different in the occultation case. When $>>1, 
we find 

$/a 1/a 1/3 4/3 -7/3 

P - / dk f dk • k - $ a . (7) 

'o X } o y 

Thus the effect of flattening either the stellar disk or the Fresnel zone by o 

4/3 

factor 4'>>1 is to increase the mean-square scintillation by a factor of 4* . 

How let us return to Eq. (1) and put in the integral over dz. Since the whole 
planetary atmosphere is at a distance z «* D that is large compared to the planetary 
radius, we can approximate 

/dz • C n 2 (z)z 2 i .z, ...) s D 2 F(D, ...)/c n 2 (z)dz, (8) 

where /c 2 dz, the integrated turbulence along the line of sight, behaves like an 
air-mass factor. For example, for an atmosphere of radius R with an exponential 
distribution of turbulence with scale height H , we have (Young, 1969b; Fabry, 1929) 

/c n 2 (z) dz - C n 2 (0) (2nR H )** exp (-2h/H fe ) , (9) 

2 

where C n (0) is the turbulent structure parameter at some reference level in the 

atmosphere where h *» 0, and li is the minimum height of the ray above this level. 

2 

For an isolated turbulent layer of thickness b and strength , such as Woo and 

Ishinnru (1974) have assumed, we have 

fc 2 (z)dz - 20 2 (2bll + b 2 ) ?i (10) 

J n n * 

= 2C 2 (2bR) ls , 
n ’ 

if the lowest point on the ray touches the bottom of the layer. The first part of 
Eq. (10) (which was first derived by Lambert in 1760) is equivalent to equations 
2 and 3 of Woo et al« (1974). If the ray passes a depth d below the base of the 




I 
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2 h 

which reduces to Eq. (10) as d *► o, and becomes b(2R/d) when d > b. Notice 
that the effective path length in the atmosphere is always proportional to the 
geometric mean of a characteristic vertical dimension and the radius of curvature. 

If we call this length L(h), we can summarize Equations (9-11) as 

Jc n 2 U)<lz ' C n 2 nax L(h), (12) 

2 2 2 
where C the peak value of C along the ray, is C (0) exp (-2h/H # .) for the 

n niaA n n l 

exponential model, and 

L(h) ** (2»RII t ) If < 13a ) 

for the exponential model, and 

L(h) = [4 R(b - h))’* , h > o , (13b) 

L(h) = [4 Kb] 11 , o ? h ^ - b/2 , (13c) 

and 

L(h) = [-2R b 2 /h] ! ‘ , h < -b/2 , (13d) 

for the one-layer model. (The approximations in Equations (13c, d) slightly over- 
estimate L(h) j at h - -b/2, the error is on the order of 27%. ) 

Now let us combine these results. Tor simplicity, assume a simple exponential 
atmosphere (Pannekoek, 1903; Fabry, 1929) with refractivity scale height 5I V . In 
the absence of turbulence, the normalized light curve <* (t) is related to the 
minimum ray height h by 

$ = <(j * “ 1 + exp (-h/H y ) (14) 

if h *» 0 where $ = this is the conventional "occultation level." We remark that 
h($) - -h (l-$); thus, for example, the 10% and 90% intensity points are 2.2 scale heigh 
below and above the occultation level, the 5% and 95% points are 2.9 scale heights 
from this level, and the 1% and 99% points are 4.6 scale heights from it. Since 
typical photometric errors are a few per cent, we obviously learn nothing reliable 
about any part of the atmosphere more than a few scale heights from this level. 


It is also easy to show that the vertical component of ray velocity through 
the atmosphere is directly proportional to the light remaining: 

- v y (*) « $ v v (l) . (15) 

If we measure time from the half-intensity point on the light curve, v y (t) -► H^/t, 
for large t. These simple relations provide some ’’feel" for the occultation 
phenomena. 

Woo and Ishiwaru (1974) have assumed that v y is the only important component 
of the motion. This is true only for central occupations. In general, there is 
also a horizontal component parallel to the planet's limb. As long as the 
occultation point does not move far along the limb, we can regard as constant; 
this is generally satisfactory for stellar occupations. For radio occupations 
of spacecraft, a large motion along the limb occurs. As Link (1969) points out, 
the motion along the limb is the projection of the geometric position of the source 
in the shy from the center of the planet's disk to the limb. Thus can be 

strongly accelerated, especially near the middle of a nearly-ceutral occultation. 
Because of the rapid decrease in v and the increase, in v^, the motion is pri- 
marily horizontal for the deeper parts of every occupation. Thus, the assumption 
of Woo and Ishitnaru (1974) is strongly violated by the data (Woo et al. , 1974) 
to which they have attempted to apply their theory. 

As long as the scintillations are weak, we can treat the observed light 
curve as an ideal one modulated by the scintillations. The modulation power is, 
from Equations 7-13, 

**/3 _7/3 2 2 

P - * a C n (0) L (h) LT , (16) 

where the major axis a of the averaging regior. is 

O dlff = ( m h (17a) 

for Fresnel-zonc averaging, and 

a * “ rD (17b) 

if the stellar disk is larger. If they are comparable, they can be added 
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quadratically (Young, 1969). Note th »t 4 depends on h and the refractivity scale 
height ll v through Eq. (14), and that L depends on h and the vertical structure of 
turbulence through Eq. (13). 

We are now in a position to estimate where, in various occultaticn curve*, 
scintillation should appear. Wc first remark that P must generally incre^w 
rapidly with depth in the occultation, not only because of the rapid increase of 
but also because of the expected increase of J C^**dz with depth. Even in 
the single-layer model of turbulence, where L decreases slowly (with about the 
square root of ray depth below the layer) , the exponential increase in 4' with depth 
must overwhelm the small decrease in L. Thus, the reason offered by Woo et al. 
(1974) for using a layer model of turbulence cannot be correct. This error is a 
direct result of omitting the dispersive flattening of the Fresnel zone *n the 
theory by Woo and Ishituaru (1974) . 

He J -_ t us estimate P, assuming that C is proportional to density, and 
adopting the following values for stellar scintillations produced in the zenith 
by the Earth's utmospher* (Young, 1967, 1969, 1970b) : P ~ 0.2 for a c 5 cm ** 

5 x 10~^ km, 4 s 1 (no refractive flattening), L { D : 8 km, number density 
p ~ 10 woleeules/ce. For a planet, we shall assume an exponential atmosphere 
with 11 v - 11 1 *» 11, adopt Eq. (13a) for L(h), and call che number density Op. We 

I 

can ignore the effects of small differences in temperature and molecular refrac- 
tivity between planets; these are generally only a factor of 2 or so. We then 
have 



for pp measured at the minimum ray height h, and all dimensions in kilometers. 

The circumstances of several occupations, and the corresponding values of P, are 
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given in Table XI for <j> « 0.5, 0.2, 0.05, and 0.01 (4*2, 5, 20, and 100, 

respectively.) The relative values of P at these levels, just due to the varia- 

4/3 2 46 

t'on of $ pp (h) with height, are 2.51, 137, 2x10 , and 4.5x10 , respectively. 

The very rapid increase in scintillation with Jcpth is due to the combined effects 

of increasing atmospheric density and Fresnel-zonc flattening} the amplitude of the 

fluctuations increases by a factor of over 1000 between the 4> *» 0.5 and $ « 0.01 

levels. 

A few special circumstances should be noted in Table II. The star 6 ScoA is 
double; this decreases P by about a factor of 2, which I neglect in the present 
crude treatmen . The appropriate value of D for a spacecraft occultation is the 
spacccraf t-to-planet distance instead f the Earth- to-planet distance, as pointed 
out by Golitsyn and Tatarski (1972) . We may also note that a_ is on the order of 

3 

a kilometer, and L on the order of 10 km, *or all the events in Table II; the 

factor a ^ J L is similar for all. Thus, because the attenuation ■'* is proper- 

2 2 

tional to ppD, the value of ppD (or pp D", which appears in the scintillation 
equation 18) at = 0.5 is similar for all planets. Thus, the value of P Is not 
greatly different from one planet to another at the same point in the light curve 
($ = constant). 

Scintillations are first visible above the noise when their amplitude exceeds 

-3 

a few per cent, say V - 10 . This can be expected near, or shortly after, the 

reference level $ = 0.5, if other planets’ upper atmospheres are as turbulent as 

— 2 

ours. When Vi? 10 > the RMS scintillation exceeds 10S, and an intensity record 

shows asymmetrical excursions ("spikes’ 1 ) because the logarithm of the intensity 
(rather than the intensity itself) has a normal distribution (Tatarski, 1961 f 
1967; Ochs and Lawrence 1969; Young, 1969, 1971); this occurs by the time the 
average intensity has fallen by a factor of about ten. On the extreme tail of 
the light curve, the theory predicts P > 1; I shall say more about this in 
subsection 4, below. 
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2, Temporal Power Spectrum 

Now let us consider the temporal power spectrum of the fluctuations. As 
long as » e scintillations are "weak" (i.e., P < 0.1 or so) this can be calculated 
by well-known methods (Reiger, 1963; Young, 1967, 1969; Lee and Harp, 1969), out- 
lined above. But now we must allow for both the flattening of the Fresnel zone 
and the flattening of its trajectory through the atmosphere (see Fig. 1) . 

For the moment we regard the planetary limb as straight. Since the cutoff 
frequency in the temporal power spectrum is 


and the cutofi. wavenumber k “ 

c 


f 

c 


V, k 

c 

2it 


1/a when 4> K 1, 


J 

we have 


(19) 


f 

c 


« V 

o 

2ita 


( 20 ) 


at the top of the atmosphere, where V is the initial speed of the ray relative to 
the planet. The horizontal components of Vj^ and £ are unaffected by differential 
refraction; the vertical compression of _a to _a4> is compensated by the reduction of 
the vertical ray velocity by the same factor, when the light is reduced to a frac- 
tion $ of its original intensity. Hence, to a first approximation, the scintil- 
lation frequency spectrum remains grossly similar throughout the occultation. 

However, as can be seen in Fig. 2, the effect of flattening the Fresnel zone 
is to stretch the region of k-space that contributes substantial scintillation 
into an ellipse. The spatial frequency components that contribute to a temporal I 
frequency f now’ lie along a line that cuts the ellipse obliquely. Thus the de- j 

jl 

tailed shape of the frequency spectrum is altered; it is no longer given by the j 
expressions derived by Woo and Ishimaru (1974) , even apart from their omission of j 

/ / n j 

the factor t> f ' . Qualitatively, the effect is to smear out the fine structure 1 
(oscillations) in the curve, and to make the knee at f c less pronounced and more ! 
gradual. ] 


If we consider the deeper phases of on occultation, where limb curvature 
cannot be neglected, we see that the effects of projection enlarge the horizontal 
extent of the Fresnel zone in proportion to the increased horizontal velocity 
component. Thus the cutoff frequency, again, is unaltered, 

A special remark must be made about Vj_, however. Woo et al. (1974) have 
allowed this velocity to be a free parameter, and have adopted values of a few 
hundred meters per second. But, as pointed out above, the value V Q (which is 
fixed by orbital motion) must be used in Eq. 20. The appropriate values of V q 
and f appear at the bottom of Table 11 for several events. For the Venus/M5 
occultation, V o is about 20 times larger than Woo et al. (1974) assumed. In fact, 
the theoretical value given in Table II (0.9Hz) is close to the spectrum turnover 
shown in Fig. 3(c) of Woo et al. (1974) for this event. 

Below this spectral knee, the power density is on the order of 0.1 Hz 1 for 
the Mariner 5 scintillations in a region where the average power received is about 
13 db below the unocculted signal (i.e., * * 0.05). Thus, the observed value of 
P is a little less than 0.1, while Table IX predicts V * 0.25 for this region. 

The agreement is fair; however, the slight discrepancy, together with the higher 
refractivity of CO,, compared to air, suggests that the stratosphere of Venus is 
a little less turbulent than our own. 

For the Jupiter/fi Seo event, the frequency cutoff implies a characteristic 
time scale (-a/V of 0.05 seconds, which should correspond to the half-width 
(at half-maximum) of the narrowest "spikes" in the light curve. Inspection of 
the spike profiles published by Liller et al. (1974) confirms that this is the 
case. Their light curves also show that the scintillation, became strong (i.e., 

P r 1) for 4> between 0.2 and 0.5, again in agreement with the prediction in 
Table II. 

The published records of the i!eptune/BD-17°43SS occultation are noiser. 

The scintillations have a large amplitude for * < 0.5, and their width appears 


r 
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to be on the order of 1/3 second, ns well as can be estimated from recordings 
with a rather compressed time scale. The predicted half width is exactly 1/3 
second, which is satisfactory. However, the scintillations are obviously much 
stronger than expected, which suggests that Neptune's upper atmosphere is con- 
siderably more turbulent than the Earth's. 

On the whole, the rough calculations of Table II are in good agreement with 
the observations, both the times of onset of "spikes" and, especially, their 
widths, are what one would expect from scintillation theory. Note that Table 11 
merely assumes Earth-like turbulence, with my free parameters to fit the data. 

3. Turbulence _vs. Layers 

The "spikes" have been attributed to layered structure in the occulting 
atmosphere by a number of investigators. Obviously, if one assumes the atmos- 
phere is spherically symmetrical, as the numerical-inversion people do, a spike 
in the light curve must produce a kink in the deduced thermal profile, but the 
required thermal excursions In the "layers" are large — 10° K or more. 

How can random turbulent fluctuations produce such large fluctuations in 

the light curves, when turbulent temperature fluctuations are typically only a 

few tenths of a degree (KMS) at separations on the order of 1, km? One might 

expect that the random effects of turbulence would be much less effective than 

the organized effects of a uniform layer; But the layer can only act. in the 

vertical direction, while turbulence can "focus" in two dimensions; and the ray 

path through a thin layer is much shorter than the effective path length 

I, (~ 10 km) through the turbulent atmosphere. Furthermore, a layer retracts 

a ray only twice (on entry and exit), but turbulence produces many random 

3 

deviations along the path. For L ~ 10 km and a ~ 1 km, isotropic turbulence 
produces about a thousand random refractions along the ray path. The net 
deviation is on the order of the square root of this — say 30 — times the 
deviation produced by a layer. For the same net deviation, the temperature or 
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refractive-index fluctuations can thus be ' 30 times smaller If turbulence is the 
cause, rather than layers. Thus, the roughly hundred-fold smaller temperature 
variations needed by the scintillation model are readily accounted for. 

There is a further problem with "layers". The fine structure in the light 
curves always has a time scale of n/Vj.. This is naturally explained by the 
spatial-filtering effect, if scintillation is the explanation; smaller-scale 
structure is smoothed out. But if layers are the explanation, we have two unsatis 
factory alternatives: either the layer thickness just happens to match the 

shrinking Fresnel sire at every depth, (an implausible coincidence); or there is 
still smaller detail in the layers, which is smoothed out by Fresnel filtering. 

But in the latter case, the finer structure must have still larger temperature 
excursions and gradients. However, the inferred excursions in the "layers" are 
already tens of degrees, and the gradients in some cases are close to the 
adiabatic value. The finer structure must exceed these values. If we do not 
accept temperature fluctuations of 100° or more, nor superadiabatlc g-adients, 
we are forced to accept the implausible coincidence cited above. As turbulence 
can produce the observed "spikes" with much smaller temperature excursions, the 
scintillation model avoids these uncomfortable problems. 

Freeman and Lynga (1970) and Rages et al. (1974) have argued that the good 
correlation between "spikes" seen at two telescopes 200 in apart at Mt. Stromlo 
is evidence for layers extending at least this far along Neptune's limb. But 
because the projected diameter of BD-17°43BS at Neptune was about 10 km, the 
atmospheric columns sampled by the two telescopes overlap almost completely; one 
must expect better than 95£ correlation between the scintillations observed at 
two points so close together. Hence their argument is worthless; the tiny diff- 
erence between the nearly perfect correlation expected from turbulence, and the 
perfect correlation expected from layers, is buried in the noise. Stations 
separated by at least the length _a along the limb must be used to distinguish 
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layers from turbulence. 

Rages ct al. conclude that no significant correlation was observed between 
two stations in Japan separated by 260 km, however. But this is an order of 
magnitude smaller than L, which requires the horizontal extent of the assumed 
“layers" to be longer in the line of sight than across it. Such accusing "fingers" 
pointing at the observer are generally a sign of some error. They ai*e removed by 
assuming isotropic turbulence instead of a layered structure. 

Another way to obtain two nearby ray paths is to use two different wavelengths 
(Young, 1969; Lillei* et al, 197*1). As the two l-resnel zones become more flattened 
and their vertical separation increases, deeper in the atmosphere, the two 
regions sampled at the same time eventually cease to overlap. However, if the 
motion is nearly vertical, different wavelengths will pass through nearly the same 
air column in succession, and all fluctuations will be reproduced at different 
wavelengths after some time delay. This is seen as chromatic scintillation of 
stars near the terrestrial horizon; when viewed through a prism, the scintillations 
appear to run through the spectrum (Respighi, 1372; Zwieky, 1950). In the Jupiter/ 
0 Sco occulta Lion, it was seen as a delay in the appearance of "spikes" at differ- 
ent wavelengths (Lille r et al. , 197a) . 

Now, if the path of the star outside the atmosphere makes a 45° angle with 
the planet's limb, a significant drop in the peak time-lagged correlation between 
spikes observed at two wavelengths should be found when the dispersion separates 
the Fresnel zones by their (flattened) thickness, because they are then separated 
horizontally by their diameter as they pass a fixed depth (see Fig. 3). This 
occurs when the spikes in the two light curves become separated by their own 
widths, as at B in Fig. 3. The more obliquely the star meets the limb, the more 
rapidly the maximum correlation should decrease with depth, if turbulence is the 
cause of the spikes. On the other hand, if horizontal layers are involved, the 
pattern of spikes should repeat equally well in the two light curves, whether the 


m 20 •* 


spikes are separated by large or small time delays at the two wavelengths. The 
data of Liller et al. (1974) obviously agree better with turbulence than with 
layers: B Sco A met the limb at roughly a 30° angle, and the spikes reproduce 

more and more poorly as their separation in time increases. Indeed, when they are 
well separated in time, one light curve can show spikes that do not appear at all 
in another, as is obvious from Figs. 2 and 4, and Tables II, III, and IV, of Liller 
et al. (1974). Such observations of a still-more-grazing occultation would show 
this effect still more strongly. 

However, there is even stronger evidence against layers in Figs. 10 and 13 
of Veverka et al. (1974a). Fig. 10 shows the temperature vs. refraetivity profiles 
for three different wavelengths; the profiles show similar wiggles (due to spike 
structure), but at different refraetivity levels. In Fig. 13 the same data are 
plotted as refraetivity scale height vs. depth, but have been arbitrarily "adjusted 
to bring the fine structure of the profiles into alignment." However, this adjust- 
ment is not legitimate. To see why, one needs only recall that the refraetivity of 
a transparent gas can change by only a few per cent, at most, over the spectral 
range (3330 - 0201 S) covered by these data. Thus, constant refraetivity is (within 
a few per cent of a scale height, or about 2 km on Jupiter) constant height. Hence 
the mis- match of structure inferred from different wavelengths is real, and cannot 
be removed by an arbitrary "adjustment". 

How docs the scintillation model handle the same data? The maxi mm overlap 
between the two Fresnel cones shown in Fig. 3, and hence the maximum correlation 
between "spikes" on the corresponding light curves, occurs when the bluer cone has 
reached a point between B and B' , if a red "spike" occurred at B. Hence the same 
atmospheric structure is traversed by the bluer Fresnel zone at a time earlier than 
would be expected for the layered model, i.e., before it reaches the depth h^« 


Thus the numerical inversion of these spikes will produce a feature in the "blue" 
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refractivity profile that occurs at a greater height in the atmosphere than the 
corresponding feature on the "red" profile. This is exactly what Fig. 10 of 
Veverka et al. (1974a) shows. Furthermore, if the profiles are "adjusted" to 
agree at one depth, the same discrepancy will occur, increasing with depth, at 
lower parts of the two prof^es. This is clearly evident in Fig, 13 of Veverka 
et al. (1974a). These effects by themselves are enough to reject the "layer" 
model in favor of turbulence. 

Grazing occupations offer another way to distinguish between layers and tur- 
bulence. If the star meets the limb at an angle 0, the time for the Fresnel zone 
to traverse its own vertical dimension is proportional to cosec 0, and the time 
scale for "spikes" due to layers is increased by this same factor. For the Jupiter/ 
B Sco event, the factor was about 2.1 for the data of Veverka et al. (1974a) and 
Liller et al. (1974). Thus, if layers caused the "spikes", they should have been 
twice as wide as observed; the narrowness of the spikes roc;ui res turbulent fluctu- 
ations of length a along the limb. This can be seen also from a comparison of the 
observed and computed light curves shown in Fig. 9 of Veverka et al. (197<’<a): the 

calculated spikes are systematically broader and less extreme than those observed. 

Fig. 1 of Veverka et al. (197Mb) shows the same defect, although the Neptune 
event was nearly central. This is perhaps the most telling point* of all: layers 

arc unable to explain the details of the light curves, which are readily explained 
by isotropic turbulence. 

111. Saturation of Scintillation 

When the value of V calculated from weak-scintillation theory is comparable 
to, or larger than, unity, the theory breaks down. Observational ly, the scintil- 
lation modulation "saturates" at a maximum value; beyond this, if the instrumental 
frequency response is limited, the measured signal variance may decrease (an 
effect sometimes called "supersaturation.") These, effects have been observed not 
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only with terrestrial light sources, but also in interplanetary and interstellar 
scintillations at radio frequencies, as well as ionospheric scintillation. The 
literature in each of these areas is quite large and will not be cited here; the 
phenomena are familiar to all scintillation workers, under such diverse names as 
"strong scattering", "multiple scattering", "focusing", "multipath propagation", 
etc. 

The theoretical treatment of saturation has been difficult. In weak scintil- 
lation, we can regard the intensity fluctuations as due to weak focusing and 
defocusing effects, which cause the first Fresnel zone (or the solid angle sub- 
tended by a finite source) to increase and decrease in area. In strong scintil- 
lation, the Fresnel zone breaks up into a "seeing disk" (Young, 1969, 3970) con- 

1 ,. 

taining "speckles". Thus, in place of a *= (XD) two length scales appear in the 
shadow patterns one corresponding to the overall blur (which imposes a low— pass 
spatial filter) and one corresponding to the speckles, which can be regarded as 
multiple-beam, high-order interference fringes. In broadband observations of 
stellar scintillation, both bandwidth (temporal coherence) and atmospheric- 
dispersion effects smear out the fine structure, and the overall "seeing" blur 
can be treated as a low-pass filter (Young, 1969, 3970). However, in the present 
case, the liighex'-order fringes or "speckles" must be Included in the theory 
(Clifford et al. , 1979; Yura, 1979; Clifford and Yura, 1979). 

The effect of saturation on the frequency spectrum (Yura, 1979b) ;s twofold; 
the "knee" in the spectrum becomes more gradual and shifts to lower frequencies; 
and the "tail" flattens, extends to much higher frequencies, and contains an 
increasing fraction of the total modulation power. The appearance of an intensity 
record changes from a dense population of oscillations (the band-limited noise 
character of weak scintillation) to a sparse population of isolated large spikes, 
superimposed on a more slowly varying background. This progressive change in 
character is well illustrated by Fig. 1 of Gracheva (1967); the stronger examples 
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look exactly like the tails of occultation curves. In particular, the isolated 
"spikes" in the far tails of the curves are due to the saturation predicted in 
Table II. 

The effect of saturation on the frequency spectrum is clearly seen in Fig. 5 

of Woo et al. (1974), where the "knee" has been blunted and shifted to roughly 

a quarter of the weak-scin filiation cutoff. Comparison with fig. 2 of Yura 

2 

(1974b) shows that this corresponds to a value of (or P, in our notation) of 
about 100. As region C of Woo et nl.(1974) is some 20 db down from the unoccult- 
ed signal, we have $ ** 0.01, at which Table II predicts P = GO. The agreement is 
within the uncertainties; the shape of the spectrum deep in the tail of the event 
seems to be adequately accounted for by saturation effects. 

There is thus ma evidence In the Mariner scintillation spectra for effects 
of a finite outer scale. There is also no evidence for "layers" of turbulence, 
if saturation (and "supersaturation") effects are taken into account; a simple 
model with roughly constant temperature structure function throughout explains 
all the data. 

IV. Effects of Saturation on Occultation Curves 

As saturated scintillation involves ray crossings, deeply saturated scintil- 
lations completely vitiate the numerical inversions of occultation curves, quite 
apart from the fact that scintillations produce spurious detail in the resulting 
profiles. Indeed, some of the numerical-! aver tors have already suspected ray- 
crossings in the deep parts of curves, where there are saturated scintillations. 
For example, Kliorc et al. (1974) mention that "Discontinuities in the observed 
Doppler shift indicate the possible occurence of multipath propagati n," and 
Veverka et al. (1974a) admit that "we suspect that some ray crossing occurred." 
Veverka et al. (1974b) and Vasserman and Veverka (1973) also admit that "in 
seeking a solution by inversion it is implicitly assumed that ray crossing is 
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negligible." But they then beg the question, saying, "at any rate ... there is 
no way in which ray crossing can be dealt with ... the assumption that it is 
negligible is unavoidable." 

On the contrary, it is possible to estimate the statistical effects of 

scintillations. Weak scintillation can be considered noise, and can be filtered 

out by simple least-squares model- fitting, or by explicit filtering (Kovalevsky 

and Link, 1969). Strong scintillations are log-normally distributed, which 

presents a statistical problem: least-squares fits are efficient only for normal 

(Gaussian) noise. Strong scintillations can be rendered normal by using log $ 

instead of but the problem then is that the logarithm of the mean intensity 

is not the mean logarithm of the intensity. Owing to conservation of energy, 

the scintillations cannot add or subtract energy from the shadow pattern, but 

* 

merely move it around. It is readily shown that 

(log i) - log (l) -o 2 /2 , (21) 

2 

where o - log (14V) is the variance of log I, the angle brackets denote averages, 
and all logarithms are Napierian . Thus, if the course of the scintillation strength 
P is known, the systematic correction (21) can be applied to the course of log I, 
and a model-fitting by least squares will be adequate. 

In the case of saturated scintillations, the problem is more complicated. As 
mentioned above, scintillation I'edistributes the light on the receiving plane. As 
long as the scintillations arc weak, the lateral scale of this redistribution is 
just the Fresnel zone size. But in saturation, it is the projected seeing-disk 
size, which can be much larger, and increases exponentially with occultation 
depth. On the sloping occultation curve, the effect of such a redistribution is 

*Thus, the procedure of rejecting "spikes", as done by Freeman and Lynjja 
(1969, 1970) and suggested by Hubbard et al. (1972) as "justified if the flashes 
were due to random lumps in the ... atmosphere", is an unacceptable violation of 


energy conservation. 
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to make the tail of the curve brighter (and hence, longer) than it should be. 

This corresponds to a fictitiously high scale height or atmospheric temperature. 

1 have remarked above on the anomalously high temperature derived from the 
Vcnus/Regulus and Jupiter /Pioneer 10 occulta tions. It is also well known that 
the Mariner 5 temperature profile shows higher temperatures than the direct Venera 
probe measurements, in the region below about 55 Pan altitude. This is just the 
region of super-saturation - i.e. , the onset of very strong angular scattering - 
according to the present interpretation (regions B and C of Woo ct al. , 1974) , 

Thus the Venera profiles are probably more reliable in this region. 

It may be argued that the agreement between Mariner ingress and ogress pro- 
files shows that the structure derived is correct. Actually, the agreement 
merely shows that the atmospheric structure (including turbulence profiles, and 
th~* onset of saturation) is similar for both ingress and egress. It does not 
prove that the structure derived by ignoring ray-crossings is correct. 

It may be objected that the radio profiles were derived from phase rathci 
than amplitude data. Multiple scattering makes the optical path between trans- 
mitter and receiver longer than it would be in a non- turbulent atmosphere with 
the same temperature profile; the "image" of the transmitter thus appears to be 
too far from •* *,e receiver, as if the atmosphere were less dense (warmer). 
Furthermore, the phase information is affected by the large-scale turbulent com- 
ponents, much larger than the Fresnel cone (Lee and Harp, 1969); large phase 
perturbations can occur before lat*ge scintillations. 

Finally, the kind of fluctuations in ionospheric electron density that 
cause ionospheric scintillations on Earth can also produce, ionospheric scintil- 
lations at other planets at radio wavelengths. This is particularly likely to 
be a problem at Jupiter, with an extensive and active ionosphere. Indeed, the 
complex structure derived in the lower Jovian ionosphere from Pioneer radio 
occultations may be analogous to the spurious detail derived from the £» Seo 
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occultation. A comparison with the much weaker ionosphere of Venus is instructive 
Figs. 9 - 18 of Fjcldbo and Eshleman (1969) show that moderate scintillations were 
observed at 423.3 MHz (70 cm) and saturated ionospheric scintillations were seen 
nt 49.8 till z ( 6 meters). They add that "Because of the multipath propagation at 
49.8 MHz, inversion of the Doppler data does not yield an accurate profile in the 
lower ionosphere." This remark applies equally well to the lower parts of tha 
S-band occupations discussed above. 

In short, the numerical inversion technique yields spurious results in the 
parts of occupations where strong scintillations (as indicated by "spikes" or 
multipath propagation) occur, i.e. , the lower parts of all occupations. As 
Washerman and Veverka (1973) say, "in the case of occultation curves containing 
spikes, this method cant Jt bo applied meaningfully, and attempts to do so lead to 
deceptive results. There! ore, this approach must be avoided in such cases." 
(Curiously enough, they were referring to the model-fitting method; but their 
remarks apply equally well to the numerical-inversion method.) 

X have just shown that scintillation effects vitiate the results of numerical 
inversion iti the lower parts of Decollations. But Vasserman and Veverka (1973) 
also showed this method was unreliable for the upper 4 scale heights of the inver- 
sion. As it is not practical to start the process until the intensity has fallen 
by a few per cent, their supposedly reliable part Joes not begin until below 
$ * 0.5. But this is just the region in which scintillation makes the results 
unreliable*. Thus the inversion method is unreliable everywhere. 

Another way of seeing this is to regard the scintillations as noise on the 
curve. Because isotropic turbulence in the planet's atmosphere can deflect 
light sideways as well as vertically, there is no guarantee that the fluctuations 
observed are a conservative redistribution of the energy in an ideal light curve 
without scintillations; large portions of the observed curve may be systematically 


too bright or too faint. The very close similarity of ideal light curves produced 
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by atmospheres with quite different temperatures and temperature gradients, shown 
by Vasserman and Veverka (1973) in their Figs. 5-7, shows that it is hopeless 
to extract more than an approximate scale, height (or temperature) at the $ * 0.5 
level from an occultatfon curve. 

Wusserman and Veverka (1973) have claimed that a temperature grudient can bu 
extracted by inversion from an occulta Lion curve, even in the presence of noise. 
However, their simulated noise is wide— bandwidth white noise, which is like photon 
noise but unlike scintillation. The values of a, in Table IX are on the order of 
a tenth of a scale height. Thus near $ “ 0.5 one has typically only 20 independent 
scintillations per scale height; if the KMS scintillation amplitude is roughly 10", 
as the observations indicate, one can expect about 2H photometric error averaged 
over a scale height, or about twice this error in the average slope of the curve 
(and lienee, the derived scale height) around ;* - 0.5. Deeper in the oceultation, 
the number of samples per scale height increase, but the scintillations increase 
still faster, by contrast, the 3*1 white noise tested by Vasserman and Veverka has 
about 100 independent samples per scale height near ij ^ 0.5, giving a slope uncer- 
tainty an order of magnitude lower than for real oceultation curves. It is the 
spectral density of the noise, not its total over a wide bandwidth, that is 
important . 

To cum up the results of this section; (1) both the gross structure in the 
deeper parts of temperature profiles derived from oceultation curves, arid the 
fine structure throughout them, are spurious. (2) In general, it is only possible 
to derive an approximate scale height, or mean temperature, for the upper parts of 
occul tat ions* (3) Contrary to the conclusions of Vasserman and Veverka (1973), a 
(modified) least-squares model-fitting procedure is probably more accurate than 
numerical inversion, which is highly susceptible to random and systematic errors 
due to scintillations. (4) In the region of disagreement, the Venera temperature 
profiles are probably more accurate than the Mariner profiles of f e Venus atmos- 
phere. 



- 28 - 


V. Discussion 

The general conclusion, that the fine structure of occultation curves tells 
us something about turbulence, but almost nothing about the general atmospheric 
structure, is a dismal one indeed. To those who have a substantial professional 
investment in the deconvolution of such details, the suggestion that they have 
wasted their time in the interpretation of random noise .cust be especially re- 
pugnant. Can we offer them hopes of either improved observational techniques in 
the future, or of obtaining more reliable information from existing occultation 
records? 

1 have already suggested a means of re-analysing existing occultation curves; 
fit aXi isothermal curve by least squares to the part of the record containing 
only weak --scintillations. This could be extended somewhat into the region of 
stronger scintillations by an iterative process: estimation of the JiMS scin- 

tillation noise (deviation from the initial isothermal curve) on short sections - 

say, each half a scale height in extent - followed by estimation of the course of 
2 

P or tt , which would be used to find the systematic corrections to the course of 
lug 1, and an improved least-squares fiL. In some favorable cases, it may be 
possible to determine a mean temperature gradient in this way. 

In the future, the largest possible telescope apertures should be used to 
observe stellar occultations, to reduce photon noise and terrestrial scintillations. 
If optical arrays several hundred meters across could be used, some aperture aver- 
aging of the planetary atmospheric shadow pattern (which has characteristic dimen- 
sions of £ parallel to the planet’s limb, and Jji perpendicular to it, in the re- 
gion of weak Scintillation) may be possible. Observations -at longer wavelengths 
would also help, because P (weak) is proportional to ?« ^ ^ (Tatarskii, 1963, 1967; 
foung, 1970). A wide optical bandpass will reduce both photon noise and scin- 
tillation (through the effects of atmospheric dispersion - see Young, 1969). 
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Alternatively, one could combine observations made at several observatories, 
and solve for the single atmospheric scale height (and possibly, its gradient) 
that best reproduces all the observations simultaneously. This seems a promising 
possibilty for the Jupiter and Neptune events. 

The spacecraft problem is more severe. If one goes to lower frequencies 
than S**band , ionospheric scintillation becomes a severe problem (Fjeldbo and 
Esbleman, 1969). On the other hand, at higher frequencies the Fresnel rone is 
smaller and the scintillations in the neutral atmosphere Increase with the 7/6 
power of frequency. Furthermore, projection of the Fresnel zone back to Earth 
gives dimensions on the order of 10 1 km for the Mariner events at Venus and the 
Pioneer events at Jupiter, so all parts of the Earth see essentially the same 
scintillation, and ground-station diversity buys no averaging of even weak scin- 
tillations. Only a receiving station on a (high) satellite can help here. 

On the other hand, spacecraft scintillations are reduced as the 11/6 power 
of the flyby distance, so closer spacecraft encounters could help substantially. 
Nevertheless, when we contemplate the expected exponential increase of scintil- 
lation with depth in the atmosphere, we must expect only modest gains in depth 
penetration. 

VI. Some Loose Ends 

The whole scintillation theory used' above assumes that the outer scale of 
turbulence is large compared to the length ji, which is typically a kilometer. 

This assumption is directly justified by the agreement between predicted and 
observed "spike" widths, er frequency turnover: if the outer scale <?ere smaller, 

the spikes would be narrower than observed, and the knee in the frequency spect- 
rum would be at a higher frequency than observed. Furthermore, Woo and Tang 
(1975) have used the wavelength dependence of scintillation to show that the 
outer scale is "large"; as both wavelengths are equally affected by the factor 

£,/3 

4> , this demonstration is not invalidated by their neglect of Frosnel-zone 
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flattening* Furthermore, there is considerable direct evidence for a Kolmogorov- 
like turbulent spectrum in our own lower atmosphere up to scales of tens of 
kilometers (Reiter and Foltz, 1967). Of course, the spectrum becomes aniso- 
tropic at scales comparable to the height ahrve ground. But the heights of 
occultation regions in planetary atmospheres are on the order of 100 km above 
the surface, so anisotropy should be unimportant. 

One may also ask whether it is plausible to expect turbulence at the great 
heights of stellar occultations. In the Earth's atmosphere, turbulence is 
strong enough to overcome diffusive separation up to heights of 100 km, and tur- 
bulence has been reported as high as 110 km (Lloyd et al. , 1973). At this alti- 
tilde the air density is less than 10 of that at sea level and comparable to 
the density at the Neptune occultation level. All the other events occured at 
higher densities. Voverka ct al. (197nb) have also argued that the occultation 
level on Neptune is near the turbopause. Thus there are reasons, apart from 
the observed scintillations* to expect turbulence in the regions probed by 
occultations. 

Wu may also inquire into the assumption that turbulence is exponentially 

distributed, with a scale height IK ~ U 4 (or approximately constant C ) . Reiger 

c v ’ t 

(1963) introduced this model, and showed it agreed with both airborne turbulence 
measurements and stellar scintillation. -Later observations of planetary scin- 
tillation (Young, 1969) supported this model, as have more recent observations 
of stellar scintillation (Young, 1970b). Modern balloon measurements (Bufton 
ct al. , 1972) also show that although there is much patchy fine structure, the 
lower atmosphere follows this trend quite well. 

VII. Future Needs 

Obviously, a more accurate cal :ul alien of the effects of turbulence on occul- 
tation curves needs to be made. We need an accurate theory for the frequency 
spectra, for example. A careful investigation of the effects of multiple 


scattering on the Mariner phase data is needed, taking account of the phase - 
locking receiver characteristics. The theory of saturated scintillations should 
be applied to the tails of occultation curves; perhaps it will be possible to 
salvage some useful information from them after all. As is suggested above, the 
existing data should be re-analysed to separate the general atmospheric structure 
from the turbulence. 

If we find that we know less about planetary atmospheres than we thought we 
did, at least we have our hands full of unfinished work; which should ultimately 
lead to a better understanding of atmospheric structure than we have at present. 

VIII. Acknowledgements 

1 thank A. Kliore for first pointing out Tresnel-zone flattening to me several 
years ago, Richard Woo has spent several hours discussing both theory and obser- 
vations with me; his recent publications prodded me into writing this paper. Carl 
Sagan has provided general encouragement. 

This work was supported by NASA Contract Number NCR 4<*i-001-117 . 


- 32 - 


REFERENCES 


Baum, W. A. and Code, A. B. (1953) , A Photometric Observation of the Occulta- 
tion of a Arietls by Jupiter, Astron. J . , 58 , 108-112. 

Bufton, J. L, , Minot t, P. 0., Fitzmauricc, M. W. and Titterton, P. J. (1972), 
Measurements of Turbulence Profiles in the Troposphere, J. Opt. Soc. Am, , 

62, 1068-1070. 

Clifford, S. F. , Ochs, G. R. and Lawrence, R. S. (1974), Saturation of Op- 
tical Scintillation by Strong Turbulence, JOSA , 64 , 148-154. 

Clifford, S. F. and Yura, H. T. (1974), Equivalence of Two Theories of Strong 
Optical Scintillation, JOSA , 64, 1641-1644. 

Fabry, Ch, (1929) , Le role des Atmospheres dans les Oceultations par les 
pi one Les, J . Obs . , 12 , 1-10. 

Fjeldbo, G. and Eshleman, V. ?.. (1969), Atmosphere of Venus as Studied With 
Mariner 5 Dual Radio-Frequency Occultation Experiment, Rad i o Sci . , 879- 

897. 

Freeman, K. C. and Lynga, G. (1969), The Occultation of BD-17°43S8 by Hep- 
tune. III. Discussion, Proc. A. S. A ., _1, 203-204. 

Freeman, K. C. and Lynga, G. (1970), Data for Neptune from Occultation Ob- 
servations, Ad . J . , 160 , 767-780. 

Golitsyn, G. S. and Tatarskii, V. I. (1972), Estimates of the Intensity of 
Turbulence in the Atmospheres of Mars and Venus, in "Physics of the Moon 
and Planets", ed. by D. Ya. Martynov and V. A. Eronshten, Nauka, Moscow; 
pp. 394-397. (NASA Technical T“anslation TIF-15,431) 

Gracheva, M. E. (1967) , Investigation of the Statistical Properties of Strong 
Fluctuations in the Intensity of Light Propagated through the Atmosphere near 
the Earth, Izv. VUZ. Radiofiz. 113, 775-7S8. English Translation: Radiophys. 
Quant. Electron. 10, 424-433. 


- 33 - 


Hubbard, W. B. t Nattier, R. K. , Evans, D. S., Tull, R. G., Wells, D. C. , 

van Cittcr* G. W. , Warner, B. and Vanden Bout, P. (1972), The Oecultation 

A 

of Beta Scorpii by Jupiter and Io. X. Jupiter, Astron. J. , 77 , 41-59. 

Kliore, A., Cain, D. L. , Fjeldbo, G. , Seidel, B. E. and llasool, S. I. (1974), 
Preliminary Results on the Atmospheres of Io and Jupiter from the Pioneer 10 
S Band Oecultation Experiment, Science , 183 , 323-324. 

Kovalevsky, J. and Link, F. (1969), Diameter, Flattening and Optical Pro- 
perties of the Upper Atmosphere of Neptune as Derived from the Oecultation 
of BD-17°4388, Astron. S An trophy s . , 2_, 398-412. 

Lee, R. W. (1969), A Review of Line-of-Sight Propagation Studies of the 
Small-Scale Structure of the Atmosphere, Radio Sci. , _4, 1211-1233. 

Lee, R. W. (1974), Remote Probing Using Spatially Filtered Apertures, JOS A , 
64, 1295-1303. 

Lee, R. W. and Harp, J. C. (1969), Weak Scattering in Random Media, With 
Applications to Remote Probing, Proe. I ELL , 57 , No. 4, April 1969, p. 375. 
Link, F. (1969), "Eclipse Phenomena in Astronomy" , Chapter 4, Springer- 
Verlag New York Inc. 

Lloyd, K. II. , Low, C. H. and Vincent , R. A. (1973), Turbulence, Billcws and 
Gravity Waves in a High Shear Region of the Upper Atmosphere, Planet. Space 
Sci. , n, 653-661. 

Ochs, G. R. and Lawrence, R. S. (1969), Saturation of Laser-Beam Scintilla- 
tion under Conditions of Strong Atmospheric Turbulence, J .O.S.A. , 59, 226- 
227. 

Osawa, K. , Ichimura, K. and Sliiraizu, H. (1968), Oecultation of BD-17°43S8 
by Neptune on 7 April 1968 (11). Observation at Okagama Station and the 
Scale Height of Neptune, Tokyo Astronomical Bulletin , Second Series, 

No. 184, pp. 2183-2185. 


- 34 - 


Pannekouk, A. (19Q3) , Uber die Erschcinungen, welche bei einer Sternbedeckung 
durch einen Planeten auftreten, Astron. Nachr. , IGA , No. 3913, 6-10. 

Rages, K. , Veverka, J. , Wasserman, L. and Freeman, K. C, C 19 74) , The Upper 
Atmosphere of Neptune; An Analysis of Occultation Observations, Icarus , 23 , 
59-65. 

Rciger, S. H. (1963) , Starlight Scintillation and Atmospheric Turbulence, 

68, 395-A06. 

Reiter, E. R. and Flotz, H. P. (1967), The Prediction of Clear Air Turbu- 
lence over Mountainous Terrain, J. Appl. Met. , 5A9-556. 

Respighi, L. (1872) , Bull. Assoc. Franc. Avanc. Sci. , jl, 1A8. 

Tatarskii, V. I. (1961), "Wave Propagation in a Turbulent Medium," translated 
by K. S. Silverman, McGraw-Hill, New York. 

Tatarskii, V. X. (1967), "The Effects of the Turbulent Atmosphere on Wave 
Propagation," Nauka, Moscow; English translation published in 1971 by 
National Technical Information Service, Springfield, Va. (TT-68-50A6A) . 
Veverka, J. , Elliot, J., Sagan, C. , Wasserman, L. and Liller, W. (1972), 
Jupiter Occultation of Beta Scorpii: Are the Flashes Time-Symmetric?, 

Nature , 240, 344-345. 

Veverka, J. and Vasserman, L. (197A), The Regulus Occultation Eight Curve 
and the Real Atmosphere of Venus, 1 earns , 21 , 196-198. 

Veverka, J. , Wasserman, L. 11., Elliot, J., Sagan, C. and Liller, W. (1974), 
The Occultation of 6 Scorpii by Jupiter I. The Structure of the Jovian 
Upper Atmosphere, A. J. , 79 , 73-84. 

Veverka, J. , Vasserman, L. and Sagan, C. (1974b), On the Upper Atmosphere 
of Neptune, Ap . J . , 189 , 569-575. 


Wasserman, L. and Veverka, J. (1973) , On the Reduction of Occultation 
Curves, Icarus , 20 , 322-345. 


- 35 


Woo, R. and Xshimaru, A, (.1974) , Effects of Turbulence in a Planetary At- 
mosphere on Radio Occultation, IEEE Trans. Ant« Prop. , AP-22 . No. 4, July 
1974, pp. 566-573. 

Woo, R., Ishimaru, A. and Kendall, W. B. (1974), Observations of Small- 
Scale Turbulence in the Atmosphere of Venus by Mariner 5, J. Atmos. Sci. , 
31, 1G9S-1706. 

Woo, R. and Yang, F. (1975), Dual Frequency Observations of Turbulence in 
the Atmosphere of Venus by Mariner 10, Bull. A. A. S. , (in press) . 

Yura, H. T. (1974a) , Physical Model for Strong Optical-Amplitude Fluctua- 
tions in a Turbulent Medium, JOSA , 64 , 59-67. 

Yura, U. T. (1974b) , Temporal-Frequency Spectrum of an Optical Wove Propa- 
gating Under Saturation Conditions, JOSA , 64 , 357-359. 

Yura, H. T. (1974c), Supersaturation: Effect of a High-Frequency Cutoff on 

Strong Optical Scintillations, JOSA , 64 , 1211-1214. 

Young, A. T. (1967), Confirmation of Reiger’s Theory of Scintillation, 

A. J ♦ . 72, 747-753. 

Young, A. T. (1969a), The Temporal Power Spectrum of Scintillation, App 1 . 
Opt. , J5, 669-SS5. 

Young, A. T. (1969b) , High-Resolution Photometry of a Planetary* Atmosphere, 
Icarus , 1 1 , 1-23. 

Young, A. T. (1970a) , Aperture Filtering and Saturation of Scintillation, 
JOSA , 60, 24S-250 . 

Young, A. T ''19706) » Saturation of Scintillation, JOSA , 60 , 1495-1500. 
Young, A. T. (1971), Interpretation of Interplanetary Scintillations, 

Ap. J. , 168 , 543-562. 

Young, A. T. (1974), Refraction Effects, in "Methods of Experimental Physic 
Vol . 12, Astrophysics, Part A: Optical and Infrared" (ti. Carleton, ed.), 


p. 139, Academic Press, 


New York. 



- 36 - 


Zwicky, F. (1950) , Seeing, Pub. Astr. Soe. Pacific , C2 , 150-155. 


37 - 


TABLE X 


Notation usad for occultations by various authors. 
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TABLE XI 


Scintillation parameters of some oecultations. 


PLANET /SOURCE 
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JUP./P10 

Scale height, 
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length, L 
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(Note: All length! 

s are in kilometers.) 



V (km/sec) 
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FIGURE CAPTIOUS 


Figure 1. (a) Path, and shape of Fresnel zone, for a 45° entry angle, (b) Cor- 

responding light intensity Cor vertical velocity component) . An 
exponential atmosphere is assumed. The lie izontal velocity component 
is constant, so horizontal position is proportional to time. 

Figure 2. Effects of Fresnel-cone and trajectory flattening on spatial and tem- 
poral frequency spectra: (a) at the start of an occultation; and 

(b) at the level s*< ** 0.5. An initial entry angle of 45° is illus- 
trated. Spatial frequencies along the dashed lane all contribute to 
the s atao temporal frequency f in both cases. Note that the Fresnel 
zone and the vertical component of V, are both compressed by the same 
factor <*.. 


Figure 3. The Fresnel cones, and their paths through the atmosphere, for two 

wavelengths. Solid lines: red; dashed lines: blue. At the beginning 

of occultation (A) the two Fresnel cones overlap completely and the 
same scintillations appear at both Wavelengths simultaneously. At 
some later point (I) the two Fresnel cones partly overlap but are 
separated by atmospheric dispersion. The overlap of the blue (dashed) 
Fresnel zone with the red (solid) cone remains large until B^, when it 
reaches the saute height h^. Thus if a spike appears in the red curve 
at B, it will also appear in the blue curve at B ' ; however, the two 
spikes may have different heights, as the overlap is only partial. 

In the tail, of the curve (C) the two cones never overlap, and scin- 
tillations become uncorrelated. But a layer of height h 0 would still 
be encountered by the blue Fresnel cone at c' ; so the correlation of 
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spike heights would still be aa good as at A, though after a tine lag. Hie two 
paths should be scaled vertically in the inverse ratio of the atmospheric rs- 
frnctivities at the two wavelengths. 


-4-2 0 2 4 6 0 10 12 

TIME— r- 



-4 -2 o 2 4 6 8 10 12 

TIME -* 


FiC.l 



a) ABOVE THE ATMOSPHERE 


FRESNEL ZONE k- SPACE 



b) AT c/> =0.5 : 

FRESNEL ZONE k- SPACE 



/ 

/ 

/ 


Tie;. 2 




